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A group of N trunks serves calls arriving in a renewal process, and lost 
calls are cleared. The number, N(k), of trunks found busy by the kth arriving 
customer is studied as a Markov process imbedded in a (usually) non-Markov 
process N(t), the number of trunks busy at t. Results of C. Palm and F. 
Follaczek on the distribution of N(k) are generalized, and a study is made of 
bounds for, and approximations to, the probability of loss. The probability of 
loss is studied as a functional of the inter arrival distribution function, and 
certain extremal properties are proven. Formulas for the mean of N(k) and 
for the covariance function are given, together with equilibrium curves for 
the probability of loss, for the mean and variance of N(k), and for the first 
four values of the covariance function. Some applications to switch counting 
arc discussed. 

I. INTRODUCTION 

We shall study a mathematical model for the random behavior of the 
occupancy of trunk groups. The principal results are complete descrip- 
tions (in principle) of (a) the variations of the traffic in time, (b) the 
equilibrium probabilities and (c), the covariance function of the traffic 
found by arriving customers. These mathematical results have practical 
application in engineering trunk groups to have a given probability of 
loss, and in estimating the sampling error incurred in certain ways of 
measuring traffic. 

A "trunk group" is a set of transmission channels (trunks) between 
central offices. The trunks in a group are often equivalent in the sense 
that a call handled on one idle trunk could as well have been assigned 
another. A "holding time" of a trunk is a length of time during which it 
is continuously unavailable because it is being seized and used as a talk- 
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ing path. By "interarrival times" we mean the time intervals elapsing 
between successive epochs at which attempts are made to place a call on 
the trunk group. With these definitions in mind, the theoretical model 
we use to describe the trunk group involves four assumptions: 

i. The holding times of trunks are independent random quantities 
having a negative exponential distribution, with mean value, 7 -1 (7 is 
the hang-up rate). This means that if a trunk is in use at time x, the 
chance that it is still in use at (x + dx) is 1 - ydx - o(dx), o(dx) 
denoting a quantity of order smaller than dx, irrespective of how long 
the trunk has been in use. The probability that a holding time is less 
than t is then 1 - exp { —yt) for t ^ 0, and otherwise. 

ii. The interarrival times of calls are independent positive variates; 
each has the general distribution A(u), where A(v) is arbitrary except 
for the condition A (0) = 0. If t k and 4+i are successive arrival times, then 

Pr{fc+, - t k S u\ = A(u), 

for all k, independently. This assumption covers Poisson (or completely 
random) arrivals as a special case. In accordance with the usage in the 
literature, we call a sequence of mutually independent, identically dis- 
tributed, positive variates, a "renewal process." The interarrival times 
in our model then form a renewal process. It has been shown by Palm 1 
and noted by Feller 2 that non-Poisson renewal processes arise in their 
own right in the study of overflow traffic from a trunk group, even when 
the original offered traffic is Poisson in character. 

iii. There are A T < °o trunks in the group. 

iv. Calls which find all N trunks busy are lost, and are cleared from 
the system. 

A model like the above, but without the strong simplifying assump- 
tion of exponential holding time, was studied by Pollaczek. 3 The model 
described in (i) through (iv) above has been considered by Palm, 1 and 
also by Takacs, 4 who used a functional equation. Takacs' paper was ap- 
parently written without knowledge of the prior work of Palm and Pol- 
laczek; in a recent paper, 5 Takacs thanks R. Syski for calling his attention 
to Refs. 1 and 3. The same model has also been treated by Cohen. 6 For 
convenience and unity of exposition, some of the results of these authors 
shall be rederived here, and attributed to the appropriate author as they 
arise. 

II. SUMMARY OF RESULTS 

It is natural to use the number N(t) of calls in progress on the trunk 
group at time t as an indicator of traffic; N(t) is a random step function, 
fluctuating in unit steps from to N. 
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Unless the arrivals form a Poisson process; that is, unless 

A(u) = 1 — exp{— u/m) 

for u ^ and mi > 0, N(i) is not a Markov process. However, let t k be 
the epoch of the feth arrival, and suppose that N(t k — 0) is known. Thus, 
we know how many busy trunks were found by the Ath call. Until the 
next call arrives at 4+i , the number of calls in progress forms essentially 
a simple death process, with death rate y per head of population. The 
conditional distribution of N(l k+ i — 0), given N(t k — 0), can then be 
calculated from the known transition probabilities of the death process 
(see Feller 7 ). No additional knowledge of N(t) for t < tk is of prognostic 
relevance to N(t) for / > t k , when N(l k — 0) is known. We define 

N(k) = N(t k - 0), 

where N(k) is the number of trunks found busy by the kth. arriving call. 
The variates N(k) form a Markov chain imbedded in the non-Markov 
process N(t). This Markov chain is the basic random process considered 
in this paper. 

Let the numbers a„ , n — 1, • ■ • , N be defined by 

a n = ^ e -" yu dA(u), 

Jo 

so that a n is the Laplace-Stieltjes transform of the interarrival distribu- 
tion A(u), evaluated at the point ny, where y is the hangup rate. The 
principal theoretical result of this paper is Theorem 1 in Section IV. This 
result gives formulas for the generating functions 

for an arbitrary initial distribution of iV(0). These formulas depend only 
on the numbers a\ , • • • , a N defined previously, so the entire Markov 
process N(k) depends only on these numbers. Theorem 1 determines, in 
principle, the transition probabilities of N(k) purely in terms of a t , • ■ ■ , 
a,v , and so provides a complete description of the statistical variations 
of the traffic found by arriving customers. For N(0) = 0, the formulas 
were obtained by Pollaczek;' the formulas to be given coincide with those 
of Pollaczek in this case. 

In Section V the limiting probabilities 

p„ = lim Pr|A*(A-) = n], 
k— » 

already considered by Palm, Pollaczek and Takacs, are briefly discussed. 
The quantity p„ is the equilibrium chance that an arriving customer find 
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n trunks busy; in particular, p N is the probability of loss. It should be 
kept in mind that p„ is not the probability that, if we inspect the trunk 
group at a random moment in equilibrium, we will find n trunks busy; 
the moments of inspection must be those immediately preceding arrivals. 
In Section V, also, various moments (such as the ordinary, binomial and 
factorial) and the variance of the limit distribution {p n \ are presented. 
Curves of the probability of loss, the fraction of trunks found busy by 
an arrival and the variance of {p n \ are plotted as functions of the offered 
erlangs for three choices of the interarrival distribution A{u). 

Sections VI and VII discuss bounds for, and approximations to, the 
probability p N of loss. The results of Section VI are general; those of 
Section VII are restricted to the case of regular arrivals. Consideration of 
the unrealistic (for telephone trunking) special case of regular arrivals is 
justified (in Section VIII) by the fact that regular arrivals form a limit- 
ing best case. 

In Section VIII we treat p N as a functional of the interarrival distribu- 
tion A(u). The chief results can be summarized informally as follows: 

i. For a fixed mean interarrival time and a fixed hang-up rate, the 
minimum loss is achieved when arrivals are regular. 

ii. Arriving customers can, without changing either their mean ar- 
rival rate or their hang-up rate, still make the telephone company give 
them arbitrarily bad service (high loss) by a proper choice of A (u) . 

iii. The maximum number of erlangs that N trunks can carry at a 
fixed loss probability p [the maximum being over A(u) that achieve p], 
is a number depending only on N and p. 

Section IX is a brief discussion of Pr{ N(k) = N} , the chance that the 
ftth arrival suffers loss, as a function of k. The case N = 2 is described 
in detail, and curves are included for one choice of A{u). 

Finally, Section X is devoted to the mean value E{N(k) } of N(k) as a 
function of k, and to the covariance function of N(k) defined as 

R(n) = Urn E{N(k)N(k + n)} - E 2 {N(k)}. 

k-»oo 

General formulas for both E{N(k)} and R(n) are derived, together with 
a recurrence relation for the latter to facilitate computation. The chief 
practical application of the covariance function is to theoretical estimates 
of sampling error in traffic measurement. Discussions of the use of our 
results to estimate sampling error in certain possible kinds of switch 
counting are given, together with some curves of the covariance. We 
stress that our results are for a finite, not an infinite, number of trunks. 
In particular, we show that a natural exponential approximation to the 
covariance, valid for N — °o , can be several times too large for small N. 
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III. SUMMARY OF PRINCIPAL NOTATIONS AND DEFLNITIONS 

'E' is used to denote mathematical expectation 

A r = number of trunks in the group 

7 = hang-up rate = (mean holding time) -1 

A(u) = Prjinterarrival time ^ u] 

Hi = ith ordinary moment of the interarrival distri- 
bution A(n) 



a n = c-" yu dA(u),,i = 1,2, ■•-,# 

Jo 

N(l) = number of trunks busy at t 

k- = epoch of the Mh arrival 

N(k) = N(t k — 0) = number of trunks found busy by 
the /.th arrival 

p„ = lim Pr!-V(A) = n} = equilibrium probability 

k- oo 

of finding n trunks busy 
p N = equilibrium probability of loss 

b„ = y*! ( J p„i — nth binomial moment of the 

distribution [p m \ 

M M n\b„ = £n(n ~~ 1) • • • (/i - m + \)p„, = nth fac- 
torial moment of \p m \ 

m„ = S m " Vm — 'ith ordinary moment of [p m \ 

a 1 = m-> — mi = variance of \p m ] 
PJx) = 1 + (•«• - l)e"™ 

Y 

E\x Nlk) \ = Z •''"' Pr!/.th call find m trunks busy! 

m=0 

*,(*,*) =L-''£'!.r V '"| 

A- 
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&n(z) = 2 ( Wm(2) = (ra!)" 1 X factorial moment 

m=n \M / 

generating function 
fc» = S ( ) Pr{2V(0) = m) = nth binomial mo- 
ment of initial distribution 



N-i 



Difli , as , • • • , a. v , z) = 2 ( • ) (1 _ *0 * * " (1 "" eo «) °*fi •••«#* 
U™ = n (1 - oi + a A+m ) 

m=l 

t/* (M = n (1 - Ot+m + <**+*) 

B(n) = limtf {#(*)#(* + »)} - £ 2 |iV(A:)} = covar- 

A-*oo 

iance function of iV(fc) 



I 

»i=0 



Qk = E mp m Pr{JV(fc) = tf | AT(0) = m) 



IV. DERIVATION OF GENERATING FUNCTIONS 

The behavior of a trunk group with (a) independent holding times, 
(b) independent interarrivals and (c) N trunks with lost calls cleared 
has been studied by Pollaczek 3 , who derived the generating functions 

^ z Pr { ftth arrival finds n trunks busy } , 

A- 

on the condition that the first arrival found all trunks idle. 
Palm and Takacs derived the limit probabilities 

p„ = lim Pr{/rth arrival finds n trunks busy} 

k- oo 

for the case of exponential holding times, to which we are also limiting 
ourselves here. Takacs used the equilibrium equations for the same 
Markov process N(k) as we have introduced. We shall show that his 
functional equation approach can be used to generalize Pollaczek's 
results, and to obtain further formulas of practical importance in traffic 
engineering. 
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We let P u (x) = 1 + (x — ])e~ yu . Then, by the argument of Takacs, 

E\x K(k+u \N(k)\ = f [P u (x)] 1+If(k) - S "- Nik) dA(u), 
Jo 

with the 8 symbol indicating that lost calls are cleared. Hence 

E{x Nik+1, \ = CI Z Pv{N(k) = n}P u l+ "(x) 
Jo L''^^ 

+ Fr{N(k) = N\P u N (x)]dA(u). 

Let 

*(x,z) = Zz k Eix^}, 

fc&O 

Uz) = E z k Pv{N{k) = n\, n = 0,l,--- ,N. 

A-gO 

Then tp satisfies the functional equation 
<p(x,z) = Eix™ ) \ 

+ z £° L\PM, z]P„(x) ~ Mz)[Pu N+ \x) - P u N (x)}\ dA(u). 

This is a discrete time-dependent analog of Takacs' functional equation. 
To solve it, set x = 1 + w and define the functions b n by 

Mi) = E (?)*.(*), n-0,l,...,iV\ 

Note that 

b n (z)=Mz), (2) 

<p(x,z) = £ *>.(*). (3) 

n=0 

If we now equate coefficients of like powers of w in the functional equa- 
tion (1), we obtain the following recurrence for the functions b n (z): 

b„(z) = za„ h„(z) + b n -i(z) - ( '^_ J *»(«)] + t- nil, (4) 

where 

fc = E (™)Pr{i\T(0) =m}, 

Jo 
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The terms /.•„ arc the binomial moments of the distribution of N(0), and 
represent initial conditions. Since 



it 

I 



£ PrjiV(A-) = n] = I 



for each k = 0, wo find b (z) = (1 — z) '. 
The solution of the recurrence (4) is 



b„(z) = n 



za, 



" ' " Z "> (5) 

i - ZQi 






20, 



where the first term of the products is always taken to be 1 . From this 
and (2) one can determine b K (z) and hence all the yj/„(z). The complete 
result is 

Theorem I : The generating junction $„{z) of Vv\N(k) = n\, defined by 

M*) =E z k Pr\N(k) = n\ 

fc£0 

is given by the formula 

*«(*)= g(- D'^ + ^W*), 

w/ierc f/w />„(*) o/r solutions of (4). //i. particular, the generating function 
of the probabilities that the kth arrival find all N trunks busy is 



* N {z) = bAz) = Z z k Pr{iV(A0 = N) 

k£0 



k K (l -2)(1 -**i) 



. , A-,(l - z) , ••• (1 - za N -i) 

'i'o i" ~r • • ■ ~r -^ 

, s-i 2fli z oia» ■ • ■ Q.\ 

~ Z) (N\ 1 - m , i ,/AT\ (l-«n)- ••(!-««») ' 

\ 1 / aoi W **fli a 2 • • • a* 

This reduces to Pollaczek's result (Ref. 3, p. 1470) when the system 
starts empty with N{0) = 0, since /.-,, = 1 , and JV(0) = implies that 
/,-, = for i > 0. Let us set 

r /N\ v 

Dn(%1 , & ,-■• X N , Z) = £ I ■ ) (1 ~ ZXx) • • • (1 - Z.Vj).V j+i ■ ■ ■ x„z \ 
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In this notation we can write 

Dwia-i , a> , ■ ■ ■ a N , z) = cha-> ■ ■ • a N z N [denominator of \Pn{z)\. 

Lemma 1 : The functions D N {x\ , x 2 , ■ • ■ , x N , z) satisfy the recurrence 
relations 

D N+1 (x,,- , ■ ■ ■ , x k +N , z) = zx k ^ N D N (x k , ■ ■ ■ , Xic+N-i , z) 

+ (1 - zx k )D N (x k+1 , • • ■ , x h +N , z). 

Proof of this is from the formula 

f N 



'D-OK-i)" 



V. THE STATIONARY DISTRIBUTION 

In the terminology of Feller, 7 the variates N(k) form an aperiodic, 
irreducible Markov chain; hence the limits 

p„ = lim Pr{N(k) = n) 

k- 00 

exist, and can be evaluated from the generating functions yf/ n {z) by Abel's 
theorem. The result is 

Theorem 2 : The stationary distribution of N(k) is {p n }, given by 



* =§<->>'•(" :>- 



with b = 1, and 

b. -n^-(i- p. t L N .) "n l -^\ , (« 

(I 1 — 0j { m=\ \>ii —1/0 O, J 

i i ti e i a J a 2 • • • «.V 

p N = probability ot loss = 



D.\(ai , a-2 , ■ ■ • , a N , 1) 

(7) 
_ . , , N\ l-o, , , /N\ (1 - oi) • • • (1 - a lV )^' 



ffl a 



1 / ai \N / aioz ■ ■ ■ a.x 

Theorem 2, and the loss formula (7) are due to Palm 1 and Pollaczek; 3 
these results have been rederived independently by L. Takacs, H. Scarf, 
the present author — and doubtless many others. 

The quantities b n of Theorem 2 are the binomial moments of \p n \, 
defined as 



k = t ( m ) 

m=n \n ) 



Pm, 



220 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1959 

and they satisfy the recurrence 
bo = 1, 

b„ = a n b„ + 6„_i - p N I _ A , n > 0, 

which can be solved to give formulas (6) and (7). The factorial moments 
M („) are then given by 

if 
M ( „) = nlb n = X) m(m - 1) • • ■ (m - n + l)p m , 

and they satisfy the recurrence 

M m = 1, 

M M = a n [M M + M ( „-i, - np N N(N - 1) ■ • • (N - n + 2)], n £ 1. 

In Fig. 1, the probability p* of loss has been plotted as a function of 
the average offered load, a, in erlangs, for three separate choices of the 
interarrival distribution A(u), for values of N from 1 to 8. The choices 
have been intentionally made so that the crucial quantities a„ depend 
on 7 and A(u) only via the offered load, a. The choices are as follows: 

i. Poisson arrivals are represented in Fig. 1 by a dashed line. In this 
case, a n = a/ (a -\- n). 

ii. Suppose that the times between successive arrivals are uniformly 
distributed in the interval (mi — b, m + b) for < b ^ m . The mean 
interarrival time is /xi , and a simple calculation gives 



-nyu, sinh ?iyb 

(In = g j — . 

nyb 

We choose b = m ; then a„ depends only on ymi = a -1 , and 

_„/ sinh n/a 



(8) 



a n = e 



n/a 



This choice of A(u) we shall call "uniformly distributed interarrivals;" 
it is represented in Fig. 1 by alternating long and short dashes. 

iii. Regular arrivals are represented in Fig. 1 by a solid line. For 
regular arrivals, a n = e~ nla , which is the limiting form of (8) as b tends 
to zero. 

The curve for regular arrivals (a„ = e~ nla ) always falls below the curves 
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a, OFFERED TRAFF 



Fig. 1 — The probability of loss: (i) Poisson arrivals, a„ = a/(a + n), dashed 
line; (ii) uniformly distributed interarrivals, o„ = e~ n/a (sinh n/a)/(n/a), long- 
and-short dashed line; (iii) regular arrivals, o„ = e~ n,a , solid line. 



for the other two choices. This is a consequence of Theorem 9 of Section 
VIII, according to which regular arrivals form a limiting best case, for 
which pff assumes its lower bound for fixed offered traffic a. On the other 
hand, the curve for Poisson arrivals, although always above the curves 
for the other two choices in Fig. 1, is by no means the limiting worst 
case, since there is none. For Theorem 10 of Section VIII says that, for 
given e > and offered traffic a, we can always find an interarrival 
distribution A(u) for which p N > 1 — e. 

The differences in p N for the various choices of A (u) in Fig. 1 are pos- 
sibly explainable by considering the amount of mass that A (u) concen- 
trates in the neighborhood of 0. For regular arrivals there is no mass, 
so that the system always has a "breathing spell" before the next ar- 
rival. For uniformly distributed interarrivals, there is always mass in 
a neighborhood of zero, but the density at is no larger than anywhere 
else. For Poisson arrivals, however, not only is there mass in any neigh- 
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■£-« OFFERED TRAFFIC IN ERLANGS PER TRUNK 
N 



Fig. 2 — The probability of loss as a function of load per trunk for Poisson 
arrivals, a n = a/ (a + n). 



borhood of 0, but the density is a maximum at 0, so that the damaging 
short interarrivals are, in a sense, the most likely. 

From Theorem 13 and the Palm formula (7) it can be verified that, 
as a — > =0 , the curves for the different choices of A (u) must approach 
each other and 1. But for small values of a there are substantial differ- 
ences among them. For this reason, they have been replotted in the 
separate Figs. 2, 3 and 4 as functions of a/N, the offered load per trunk. 
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£-, OFFERED TRAFFIC IN ERLANG5 PER TRUNK 
N 

Fig. 3 — The probability of loss as a function of load per trunk for uniformly 
distributed interarrivals, a,, = c~ n/ " (sinh n/a) /(n/a). 

The first two ordinary moments wh and m* of \-p n \ are respectively 
given by 



in 



i = M (1) = 6, = 53 "P« = 



Oi(l - p,v) 



1-ai ' 

A' 

m a = .1/(2) + Ma) = 26 2 + bi = 53 » 2 P« 
aia-.(l — p. v ) 2a 2 -Vp.v — wii 



(1 - a x )(l - tte) 

- a, , for JV = 1. 



1 - 09 



for N > 1 , 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 09 



TT' OFFERED TRAFFIC IN ERLANGS PER TRUNK 

N 



I.I 1.2 



Fig. 4 — The probability of loss as a function of load per trunk for regular 
arrivals, o„ = e~ nla . 

The variance associated with {p n } is then 

a 2 = fl(0) = ma - mi - 2/> 2 + &, - &! 2 , 

where 7i!(/j) is the covariance function. 

Because of the bias introduced by defining N(k) to be the number of 
busy trunks found by the Ath arriving customer, it is not in general 
true that mi equals lim E{N(l)} as t tends to <», even when this limit 
exists. In Fig. 5, the ratio 
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TOi 



~ = fraction trunks found busy 

N 

_ expect ed number found busy by an arrival 
number of trunks 

is plotted as a function of offered load a for Poisson arrivals. In Figs. 6 
and 7 the same ratio is plotted for uniformly distributed interarrivals, 
and regular arrivals, respectively. 

In Figs. 8, 9 and 10, the variance a of N(k) in equilibrium is shown 
plotted against the offered load a for Poisson arrivals, uniformly dis- 
tributed interarrivals and regular arrivals, respectively. The variance is 
also the value of the covariance function R(n) for n = 0. In all cases, 
as the load a increases, the variance increases to a unique maximum, 
and then decreases to zero. 

VI. BOUNDS FOR, AND APPROXIMATIONS TO, p N FOR GIVEN Oi, ■ ■ ■ , d N 

This section is devoted to inequalities which may be useful in esti- 
mating the loss probability p N without too much computation. Since 
1 > ai > • • • > a N , we have 

1 — a n , 1 — a n +i 

a n CIn+1 

so that, from (7), we find 

toc-i s y«--«*G)CT# 

This proves: 

Theorem 3: Tlie probability p N of loss salisjies (a N ) N ^ p N S («i) • 
To obtain a sharper result, write 

Vfl ~ l = ( ai a 2 • • • a N )" 1 2 ( . j (1 - ai) • • • (1 - a;K+i ••• a H . 

Then, in view of 1 > a* > • • • > a N , 

£ (J) (i - riWr s a -^j^ s E (f) a - *>w-'. 

From this we conclude: 

Theorem 4: 77te probability p N of loss satisfies 

(1 - a, + a*)"* ^ ^ Si (1 + ai - a,)-" 

aifl2 • • • a N 



22G 
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Fig. 5 — limju» ElNi/^l/N = nn/iV as a function of offered traffic a for Poisson 
arrivals. 
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2 3 4 5 6 

a, OFFERED TRAFFIC IN ERLANGS 



Fig. 6 — lim**«. -£7{ ^(A;) J/^ = »t,/iV as a function of offered traffic a for uni- 
formly distributed interarrivals. 
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2 3 4 5 6 7 8 

3, OFFERED TRAFFIC IN ERLANGS 

Fig. 7 — lim*-,. E\N(k) }/N = vii/N as a function of offered traffic a for regular 
arrivals. 

This result suggests that, if «i — a N is sufficiently small, then the prod- 
uct aia 2 • • • a N can serve as an approximation to p N . There are cases, 
to be exemplified later, in which this is a good approximation. However, 
the next theorem shows that the product a^ • ■ ■ a N always underesti- 
mates the loss. 

Theorem 5: For N = 1, p N = a x ; for N ^ 2, p N > a t a 2 ■ •■ a N .* 
To prove this, we write p N in the notation of Lemma 1 as 



p N = 



CLid* • • • CLtt 



D N (cii , • • • , as , 1) ' 

so that it suffices to prove that £>*(<*! , • • • , a N , 1) < 1. We shall actually 
prove the stronger result that D N (a k , ■ ■ ■ , a k+ N-i , 1) < 1 for /; ^ 1. 
First we note 

D 2 (a k , a k +i , 1) = a k a k+l + 2(1 - a k )a k+l + (1 - a*)(] - a k+l ) 
= 1 — a k + a k+1 < 1. 
Now, because 1 > a t > ■ ■ ■ > (h > • • • , we find 

D.\(ak , • ■ • , flt+A-1. jj < Av(flt +ii • ' ' » Qk+N > 1) 



dicdk+i • ■ ■ a k +\-i 



a k +ia k+ 2 ■ ■ • a k +N 



* A. J. Goldstein has pointed out that Theorem 4 implies directly that p N > 
ciidi ■ ■ • fl,v for N ^ 2. 
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2 3 4 5 6 

a, OFFERED TRAFFIC IN ERLANGS 



Fig. 8 — The variance <r 2 [= R(0)] of N(k) in equilibrium for Poisson ar- 
rivals. 
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1 2 3 4 5 6 

3, OFFERED TRAFFIC IN ERLANG5 

Fig. 9 — The variance a 2 [= R(0)\ of N(k) in equilibrium for uniformly dis- 
tributed interarrivals. 

Therefore, the recurrence of Lemma 1 gives, for z = 1, 

Ds+i(ak , ■ ■ ■ , o-k+N , 1) < D N (a k+ i , ■ ■ • , a k+N ,1) < 1, 

and the result follows by induction. 

We now discuss the approximation p N ~ a^ • ■ • a N . Since 1 > a x > 
a N , two cases in which «i — a N is small are as follows: (a) a,\ is close to 
and pn is very small; (b) a N is close to 1 and p N is very high. The quan- 
tity ai — a N determines the excellence of the approximation, as meas- 
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1 2 3 4 5 6 

a, OFFERED TRAFFIC IN ERLANGS 

Fig. 10 — The variance a- [= R(0)} of N(k) in equilibrium for regular ar- 
rivals. 

ured by Theorem 4. The value of ai — a N may be estimated from below 
in terms of dj alone by the inequality ai — a K ^ ai — ai . From Theo- 
rems 4 and 5 we see that 



(1 — ai + a*)" g r = 



aia2 • ■ ■ a.v 



< 1, 



and this inequality indicates the values of Oi — a* for which p N <~^ 
aict2 • • • «at is justified. 

To put the matter more intuitively, we note that a t is the chance that 
a conversation, in progress at one arrival epoch, is still in progress at 
the next arrival epoch, i.e., 

Oi = Pr{ holding time > interarrival time } . 
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Similarly, if hi 



, h N arc N (independent) holding times, 
-Pr< 



a N = p r J min^ ) x . > interarrival timej. 



So the approximation is likely to be good at least when the chance that 
one holding time exceeds an interarrival time is not very different from 
the chance that each of N holding times exceeds an interarrival time 
(the same one for all N). As a tentative conclusion we may say that 
Vn ~ cua-i • ■ ■ a.v is good when the loss is very high or very low. 

The ratio r = aia 2 ■ ■ ■ a N /p N has been plotted as a function of the 
average offered traffic a in Figs. 11, 12 and 13 for Poisson arrivals, uni- 
formly distributed interarrivals and regular arrivals, respectively. The 
curves bear out the conclusions of the previous paragraph, that the 
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Fin. 11 — The ratio r = ifl\(i- • • ■ on)/pn as a function of traffic a for Poisson 
arrivals. 
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Fig. 12 — The ratio r = {a&i ■ ■ • a/f)/pN as a function of traffic a for uni- 
formly distributed interarrivals. 

approximation p N "*> aia 2 • ■ • a N is good for low and high traffic. Fig. 
14 shows a detail of r for very low traffic, for all cases at once. 

Lemma 2: For m, k ^ 1, ol+i + a k+m ^ a k + a k+m +i . 

Proof: the case m = 1 holds by convexity; for the same reason, 

a k + a k+ 2 ^ 2a*+i . 
Assume that the lemma holds for a given m and all k ^ 1. Then 

j ^ Ofc + Ofc+i j 

a* + 2 t at+i+m S ~ ~r a k+2 + m , 



_l „ <- ,, _i_ ttk ~ °*+2 _|_ 
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Fig. 13 — The ratio r = (aiOa • • • aN)/vN as a function of traffic a for regular 
arrivals. 



But a k + a.A- + 2 ^ 2aA- + i implies 

, Ok — Oa+2 ^ 
«A--|1 + ^ ^ «A- , 

so the lemma follows by induction. 
Theorem G: Let 

L k {N) = ft [1 - a k + a k+m ] , 

m=l 

U k w =11 [1 - a k+j + a k+ s). 

;=0 
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Fig. 14 — Detail of r = (aia 2 • • • a N )/p N for low traffic a for Poisson arrivals, 
uniformly distributed interanivals and regular arrivals. 
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Then 

Inf*"* ^ D N (a k ,■■■ , a, + *_i , 1) ^ U^'" 
and, also, the chance of loss satisfies 

ai«2 ' " CLir < __ a i a z ' ' ' a x < fl i fl 2 • • • fly 

W-» = P ' V ~ D*(ai , fl2, • • ■ , ay , 1) = L^-» 

Proof: For N = 2, we have for fc ^ 1 : 

Di{a k , flfc+i , 1) = 1 — a* + a*+i = £4 = L k . 

Now, assume that for all k ^ 1 

L, <A '- n £ ZM«* , • • • , fl fc+JV _! , 1) g f/, (A '" n . 

Then, by Lemma 1, 

a* + »L* (Ar ~ + (1 - «*) W lV_0 ^ ZWa* , • • ■ , «*+" , 1) 

:g Ofc+r^"" + (1 " a k )U k+1 {N ~ l) . 

By convexity and Lemma 2, 

T (A'-D -„. X (AT— 1) 

ivfc+1 = ^* » 

7-t (A'-D > 77 (JV-1) 

Therefore 

L fc (A ' _1) (l - a k + flA+y) ^ D N+1 {a k , ••■ , a k+N , 1) 

^ (1 — a k + fl*+w)?7fc+i 
But 

/-. . „ \ T (N-l) j (AO 

(1 — a k + a k+N )L k = L k , 

(1 — a k + flfc-fjv)L/fc + i = Uk , 

so the theorem follows by induction. 

VII. BOUNDS AND APPROXIMATIONS "WHEN ARRIVALS ARE REGULAR 

In telephony, it is unrealistic to expect regular arrivals. Nevertheless, 
the results of Section VIII indicate that regularity of arrivals represents 
in a definite sense a limiting best case, for which the loss assumes a lower 
bound. For this reason we devote some effort to approximating the loss 
Pat in this case. 

For regular arrivals the loss ips is given by 

(N\ 1 - x 1 - x_ _ 1 - X s 

X 2 '" x> ' 



o-r-lW-i 
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where x = exp{ — (1/a)} and a = offered erlangs. A simple procedure 
for obtaining an upper bound on p N is as follows: we note that, since 
x < 1, 

(P»r l = Z (J) (i - *)'a + .ou + x- 1 + x-*) ■ • • (e x- 1 ) 
^(J)a-.w.-«i-.>-?%^. 

The term on the right of the last inequality is seen to be the reciprocal 
of B[N, 1/(1 — x)], where B(c, a) is the classical Erlang B function; 
that is, 



B(c a) = &*}} = e - 

7=0 j\ 

where P(c, a) is the cumulative term 2l a n e~ a /n\ of the Poisson dis- 
tribution. This proves: 

Theorem 7: // arrivals are regular and a erlangs are offered, then Pn ^ 
B(N, rf), where B is Erlang' s function, and t\ — (1 — x)~ l = (1 — e _I ") _1 . 

From Theorem 9 of Section VIII we know that p N ^ B(N, a) ; that 
is, we overestimate the loss for regular arrivals if we pretend that arrivals 
are Poisson. Let us therefore see whether the bound of Theorem 7 is 
better than B(N, a). Let a = (1 - f)~\ so that 17 = (1 - e^)~\ Now 
f is tangent to e f_l at f = 1, i.e., at a = 00, and e is convex; hence 
i~ x ^ f , and 1 - f ^ 1 - S~\ so that 

a = (1 - f)" 1 < (1 - e M ) _1 = r, 

for finite a. Since B is monotone increasing in the offered erlangs we 
conclude that B(N, a) < B(N, 77). Thus the bound of Theorem 7 is 
nowhere as good as the overestimate B(N, a) for p N . 

However, there is a systematic way of obtaining a useful upper bound 
on p N for regular arrivals. This bound again has the functional form of 
Erlang's formula B(N, ij). However, 77, instead of being chosen equal to 
a, is chosen to correspond to a Poisson process, which gives the right 
value of 01 , exp{ — (1/a)}, and involves fewer offered erlangs r? < a. 
Now 

fe _1 a for regular arrivals at a erlangs 
1 \*?/(l + v) for Poisson arrivals at 77 erlangs. 
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So 77 erlangs will give the right value of a\ if and only if 
, = -JL- for y = e~ lla . 

1 - y 

We first show that, if 77 is defined in this way, then 77 < a. For u > 0, 
we have u -f- 1 < e u , so that for u = cT we find 

<r ,/a < o(l - e- ,,a ), 

■n = - — - — < a for y = <T 1/a . 
1 - V 

For this choice of 77, then, B{N, 77) < B(JV, a). Now, from formula (7), 
it is apparent that if the a< are replaced term by term with quantities 
a/, with a, ^ a,', the result will be ^p N • We choose 

a/ = u/(i? + i), t - 1, 2, • ■ ■ , N. 

The a/ correspond to Poisson arrivals with 77 erlangs offered. To obtain 
a bound it remains to be shown that, f or i = 2, 3, • • • , N, 



a; = e 



V + i' 

This is equivalent to 

y + i = z.y ^ ?/ 1- ', for y = e~ 1,a , 

which is seen to be true because y -\- i — iy is tangent to y l ~ l at y = 1. 
The result of replacing a, by the chosen 0/ is just B(N, 77). This proves: 
Theorem 8: // arrivals are regular and a erlangs are offered, then 
Pn ^ B(N, 77) < B(N, a), where B is Erlang's function, and 

y e^ 

V 1 - y 1 - c<- l/a> - 

This result suggests use of B(N, 77) as an approximation to p, v . Two 
numerical cases illustrate this approximation: 

i. N = 8, 8 erlangs are offered; then y = e~°' m and 77 = 0.747. We 
find V s = 0.17, B(N, v) = 0.20, B(N, a) = 0.235. 

ii. N = 5, 8 erlangs are offered; again, 77 = 0.747, and p N = 0.437, 
B(N, tj) = 0.450, B(N, a) = 0.478. 

VIII. THE LOSS AS A FUNCTIONAL OF A(u) 

For each N, and each hang-up rate 7, the loss p N can be regarded as 
a mapping from the set of distributions A(u) of positive variates to the 
interval (0, 1). We write px(A) in this section for the loss resulting from 
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the interarrival distribution A(u), and we study the loss as a functional 
of A(u). 

First, it is instructive to keep the mean interarrival time hi fixed and 
to vary the interarrival distribution A(u). Since e~" yu is convex in u, 
we find 





— ! 

e 


™ ^ e-" yu dA(u) 
Jo 


and hence 




1 - <r wwi > 1 - a n 



But e~ n "" il = a H for the case where arrivals are regular, and mi apart. 
This proves: 

Theorem 9: If y, mi are positive constants, then 



inf 'ipx(A) I udA(u) = /xi 



A \ 

is achieved for the unit step distribution 

1, u ^ Mi , 



[0, W < Ml ■ 

Thus, the probability of loss assumes a minimum, for fixed y and mi , 
when the arrivals are regularly spaced at epochs mi apart. 

We next show that, if the mean interarrival time mi and the hang-up 
rate y are kept fixed, then the probability of loss can still be made ar- 
bitrarily close to unity by a proper choice of A(u). 

Theorem 10: If y, mi o,re positive constants, then 



sup<pjv(A) / udA(u) = m> = 1. 



To prove this, let 1 > e > be given, and consider those distributions 
which have a mass (1 — p) at ya > 0, and a mass p at jji > 0. For such 
an A (u) we have 

a n = (1 - p)<T n ™ + pj- w . 
Let q = (1 — p) exp{ — iVy?/o}, so that, for each n, 

1 - g > 1 ~ «» 

a ~ On 
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Then, since a N < a N -j < • • • < d\ < 1, we find from Palm's formula 
(7) that 



,,(-o,(. + ^f. 



We can now choose p and ij so small that p N (A) ^ 1 — e, independ- 
ently of iji , which can (hen be chosen to satisfy mi = //n(l — p) + UiP- 
This proves the theorem. 
It is natural to use 

1_ = 1 

r, 00 

/ udA(u) 

as a measure of the culling rate, and to use 

r N {A ) = ?,v = fraction served times calling rate 

Ml 

as a measure of the rate of service, the rate at which calls are actually 
being completed. Suppose now that we are willing to tolerate a probabil- 
ity p of loss. Can we find an interarrival distribution A (w) which achieves 
p and for which the rate of service is a maximum for a given hang-up 
rate 7? To answer this question, define the function 

, , 4. (N\ 1 - .r, /iV\ (l-.r 1 )---(l-.r„) 

so that p N (A) = [/(01 , <h , ■ ■ ■ , ojr)]" 1 - 
Theorem 1 1 : // 7 > and < p < 1 , then 

t /a\ 1 1 w 1 y(l - p) 
BUp{rjr(A)|pjrU) = v\ = _ j - . 

w/icre .r zs </ie unique solution of the equation f(x, x , • • • , x ) = p in 
the unit interval. The supremum is achieved by the unit step distribution 
A(u) defined by 

_ fl, u ^ -7"'log.r, 
A(u) -\0 f «< -7-lng.r. (9) 

The function /(.r, .r 2 , • • • , x 1 *) is monotone, decreasing from » to 1 
in the unit interval. Since/ is continuous, and < p < 1, there exists 
a solution x of the equation /(.r, .r'\ • • • , .r v ) = p~\ Obviously, for A(u) 
denned by (9), we have p N (A) = p. Now let B(u) be any other inter- 
arrival distribution with a finite mean, so that the service rate r N (B) 



240 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1959 

exists. Suppose that B(u) achieves the probability p of loss; i.e., that 
p N (B) = p. We show that 



f 

Jo 



it dB(u) ^ — 7 ' log x. 
For, suppose the contrary and set 

y = exp< — 7 / « dB(u) > ; 
then, by Theorem 9, [%, y\ ■■■ , y")}' 1 ^ p, and 

ii flB(u) < — 7 _l log x 



i; 



implies y > x, so that \f(x, x\ ■■■ , x N )] ' < [f(y, y\ ■■■ , y")] * ^ p, 
which is impossible. This proves that 



inf | / udB(u) 



p N (B) = p) = -7 log a;, 



and also Theorem 11. Note that the supremum in Theorem 11 is a linear 
function of the hang-up rate, 7. 

Let N(t) be the number of trunks busy at time t, and let E\N{1)\ 
be its average. It is not always true that lim E\N(l) \ exists as t — > °o. 
However, if A(u) is not a lattice distribution, then 

i-00 7Mi 

where mi niay be <» . This limit is the number of erlangs carried by the 
trunk group in equilibrium (see Takacs ) . Now a lattice distribution 
can be approximated arbitrarily closely by absolutely continuous dis- 
tributions. Thus, an immediate consequence of Theorem 11 is: 
Theorem 12: 7/0 < p < 1, and x is as in Theorem 11, then 

sup lim E[N(t)\ = 1 1 ~ P , 

A t-oo — lOg X 

where the supremum is taken over A(u) such that Pn(A) = p and such 
that lim E{N(t)} as t — * °o exists. 

This theorem means, intuitively, that the maximum number of er- 
langs that N trunks can carry at a fixed loss probability p [the maxi- 
mum being over the appropriate A{u) which achieve a loss p] is a number 
depending only on N and p. 
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It may also be of interest sometimes to know what is the least proba- 
bility of loss incurred by offering a traffic of a erlangs to N trunks, with 
A(u) being varied, and 7 as well. The answer is given by: 

Theorem 13: If a > 0, then 

inf ipAA) I f yu dA(u) = aA = [f(x, x\ • • ■ , z"T\ 

where x = e~ lla , and the inf is achieved by any unit step distribution A(u) 
and 7 > such that 

A(M) -\0, W <(a7)- 1 . 
The proof is essentially that of Theorem 9, and is omitted. 

IX. Pr(JV(/i') = N) AS A FUNCTION OF & 

The time-dependent behavior of the process JV(A;) is only touched on 
here, since a complete treatment requires the detailed investigation of 
the roots of the polynomial D N (ai , a 2 , • • • , a N , z) occurring in the 
generating function \{/n(z). Such a study is still incomplete. 

Nevertheless, some hints of the rate of approach to the limit y N can 
be obtained from Theorem 1 and \f/ N (z) as they stand. For instance, if 
N(0) = 0, then 

(1 — z)~ a\ai ■ • • a N z' 



$n{z) = 



2 ( . ) (1 — ttiz) •••(!- a,zK + i 



N-i 

a N z 



From this it can be seen directly that 
Pr{iV(fc) = N\N(0) = 0)) = 

for k < N, 
a x ai •■ ■ a N for k = N, ^ 10 ) 

aiao • • • a N 1 + S («; — «jv) for k = n + 1. 

More terms may be computed from the generating function, but the 
labor involved increases rapidly. It is to be noted that (10), together 
with Theorem 5, suggests that the approach to p N is monotone; also, 
the first nonzero term is the approximating product Oi0 2 • ■ • a N discussed 
in Section VI. 
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For N = 2 trunks, it is possible to discuss ~Pr{N(k) = N | JV(0)} in 
a particularly simple way. The results are given here, together with a 
numerical illustration, for the light they shed on the time development 
of the process. From Theorem 1 we find that 



dydiZ" 



£ z Pr [N{k) = 2 | AT(0) = 0| = "^ , , 

* (1 - z)(l - za x + za 2 ) 



so that 

Pr \N(k) = 2 | JV(0) = 0\ = 



aidi 



1 — ai + a-> 




[1 - (a, - a 2 ) fc ~ 1 ] 



for k = 0, 1 
for fc > 2. 



Here p* is Oia 2 /(l — a x + a 2 ), and is approached exponentially. 
Similarly, the generating function of PrjJV(fc) = 2 | N(O) = 1} is 



didtf 



+ 



za 2 



(1 — z)(l — za\ -f- 2a 2 ) 1 — zax + 2a 2 ' 
so that 
Pr \N(k) = 2|iV(0) = 1| = 



r/2 

dyd2 r , / \fc-ll I ( \*:-l 

, 1 11 - Ui - d 2 ) J + flkUi - a 2 ) 

1 — Oi + 02 

Finally, the generating function of Pr{iV(fc) = 2 | iV(0) = 2} is 

d\diZ l , Zd 2 



for /c = 
f or k = 1 

for k > 2. 



(1 — z)(l — zai + zaO 1 — zdi + z«2 
from which we find 
Pr \N(k) = 2 | N(0) = 2J = 



+ 1, 



fllfl2 

[ 1 - ai + a 2 



[1 — (ai — a 2 )* J ] + diidi - a 2 )' 



1 


f or fe = 


a 2 


for A; = 1 


— i 


for k > 2. 
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Fig. 15 — Pr{JV(fc) = 2\N(0) = 0) for a = 1 erlang, N = 2 trunks and uni- 
formly distributed interarrivals. 

This agrees with the previous conditional probability for k ^ 1, as it 

should. 

The three conditional probabilities Pr{iV(A-) = 2 | N(0) = m) for 
m = 0,1,2 have been plotted as functions of k for uniformly distributed 
interarrivals in Figs. 15, 16 and 17, respectively. The probabilities have 
been drawn continuously, but of course the functions are only defined 
for integers k. The example chosen exhibits a very rapid approach to 
equilibrium in terms of numbers of arriving calls, since the third arriving 
cull finds essentially the equilibrium situation. 

X. THE EXPECTATION OF N(k) AND THE COVARIANCE 

The next result gives a formula for the mean value E{N(k)} in terms 
of the initial value E[N(0)}, and the probabilities Pr(JVO') = iV} for 

j ^ k - 1. 

Theorem 14: The mean value of N(k) is 



E\N(k)} = 



ai(l — a/) 
1 - a. 



+ afEimO)) - £ «i i+1 Pr {N(k - j - 1) = N\. 
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Fig. 16 — Fr[N(k) = 2\N(0) = 1) for a = 1 erlang, N = 2 trunks and uniformly 
distributed interarrivals. 
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Fig. 17 — Prf JV(A;) = 2|iV(0) = 2) for a = 1 erlang, N = 2 trunks and uniformly 
distributed interarrivals. 
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To prove this we first obtain the generating function bi(z), either by 
differentiation from (1) or directly from the recurrences (4). We find 

s rtiiwwi = w 2 ) . "jq-*)"-'-wi + *wo>i , 

and this gives Theorem 14 upon expansion in powers of z. 

We define the covariance function R(n) of the random process N(k) 
by 

fl(n) = lim.E{ #(£)#(£ + »)} - #{#(*)}• 

From Theorem 14 we can derive a formula for the covariance function 
R(n). 

Theorem 15: // ^n.w(z) is ^s(z) for the initial condition A/(0) = n, 
\ z\ < 1, |p,„! is the stationary distribution of N(k), and w, = E ™'p m 
for i = 1,2, tten 

2^ z R(n) = 2^ Pmm 



1 — za x 






and 

B (n) = R(-n) = ^ = a '" +1) + *»m. - m t 2 

1 — Oj 

- E ^P« E «i i+1 Pr W* " 3 ~ ] ) - N I ^<°) = W J- 

Before developing the results of Theorem 15 into a form useful for 
computation, we shall sketch the reasons for interest in the covariance 
function R(n). The function expresses quantitatively the cohesiveness 
of the process, the extent to which N(k + n) and N(k) are correlated. 
Besides this theoretical role, the covariance is involved in the practical 
matter of evaluating (theoretically) the sampling error in a certain kind 
of switch count (traffic measurement.) For a concrete example, suppose 
that 

S = E w*o 

i 
is used to estimate the average traffic encountered by arriving custom- 
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ere. Here n is the number of successive observations of the random 
process N(k). The variance of S is 

varjS} = E fe £ N(i)N(j)\ - nW 

= EHoov{N(i),N(j)} 

i i 

= E Z fl(»" - i) 

= «ft(0) + 2E (n-j)R(j), 

where we have assumed that the observations began in a condition of 
equilibrium. Thus var{<S} can be expressed in terms of the covariance 
function R(?i). 

The formula for R(n) can be made more useful for computation by 
turning it into a recurrence relation for successive values of a certain 
linear function of R(n). We define auxiliary quantities Qk by 

Qk = E mp m Pr{AT(*) = N | N(0) = m) (11) 

ro=0 

and note that 

R(n) + m x - — - = ai [ m 2 - - ) - 2- «i Qn-j-i • 

1 - ai \ 1 - aj >=o 

Hence also 

BCn + D + mx 2 - mifll 



1 — ai 



= a, jar (m 2 - ^-J - Q„ - £ ai i+1 Q»-,-x} (12) 

= a, (ff(n) + W - -^L. _ q h 

[ 1 - ai 

Thus, if the Qt are known, the i?(n) may be calculated by a simple re- 
cursive procedure from R(Q), which is the variance. The calculation of 
the Qk is simplified by the fact that, for small k, (a region of principal 
interest), many terms of the sum defining Q k are 0. For example, if 
^ m < N - k, the conditional probability Fr{N(-k) = N | N(Q) = 
m\ is 0, since it is not possible for the kth man to find all trunks busy 
if the 0th man found fewer than N — k busy. The first few correction 
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Table I. — Pr {N(k) = N I JV(0) = mi. 





k 
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1 

as 
as 


2 


3 


N 


Na N a N -i - (N - 


l)aw» 


(N* - N + l)ajv 3 

/ AT* — A/\ 
+ (N - 2iV 2 + 2 Jay-iay' 

N(N - 1) 
+ Na N -i 2 aN + o aN-2a N -idf/ 


JV-1 





N* N -ia N - (N - 


Da* 2 


same as above 


N-2 








a N -ia N 


ajv_i[aAr_i + (N — l)(ajv_ 2 - CTat)) 


iV-3 











dN-2 dS-l &N 



terms Q fc as defined above may be computed (by summation) for k = 
0, 1, 2, 3 from Table I, which shows Vv{N(k) = N \ N(0) = m\, valid 
for m ^ 0: 

Curves of the covariance function R(n) for n = 1,2 and 3 are plotted 
as functions of the offered traffic a for trunk group sizes N = 2, ■ • • , 8, 
as follows: in Figs. 18 through 20 for Poisson arrivals; in Figs. 21 through 
23 for uniformly distributed interarrivals; and in Figs. 24 through 20 
for regular arrivals. The curve for N = 1 is not shown in any of Figs. 
18 through 26 because, in this case, R{n) = for | n | > (see below). 

The following conclusions seem to be reasonable after examination of 
the curves: 

i. R{n) is nonnegative and monotone decreasing in \n\. 

ii. For n and traffic a fixed, the covariance R(n) for Poisson arrivals 
exceeds the covariance R(n) for both the other two interarrival distri- 
butions (uniform and fixed) we have considered. Similarly, the covari- 
ance R(n) for regular arrivals falls below the value of R(n) for both 
Poisson arrivals and uniform interarrivals. We conjecture that R(n) for 
regular arrivals is less than or equal to R(n) for any other distribution 
of interarrivals, for the same traffic. 

A particularly simple but important case arises when N = 1 ; the case 
is simple because R(n) = except for n = 0; the case is important, not 
because groups consisting of a single trunk are common (they are not), 
but because the case A r = 1 corresponds to making a measurement only 
on the first trunk of a group (of arbitrary size) in which the trunks are 
tried in a fixed order. For N = 1 it is easy to see (from Theorem 1) that 

Fr[N(k) = 1 | .V(0)) = o Ni0hl for k = 0, 
= a-i otherwise, 
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Fig. IS — The covariance value R(l) as a function of traffic a for Poisson 
arrivals. 
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Fig. 19 — The covariance value R(2) as a function of traffic a for Poisson 
arrivals. 

so that N(k) is independent of N(0) for k > 0. Thus, in this case, 
R(0) - var [N(k)} - at - d 2 , 
B(n) = 0, for n ^ 0, 



var{iS/n} = 



Oi - Oj 



so that 5/n is a consistent and unbiased estimator for a% . It is to be 
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1.8 



R(3) 




2 3 4 5 6 

a, OFFERED TRAFFIC IN ERLANGS 

Fig. 20 — The covariance value 72(3) as a function of traffic a for Poiason 
arrivals. 

emphasized that in this case S is the sum of n independent identically 
distributed random variables, each equal to 1 with probability a x , and 
to with probability 1 - a x . Thus, S has a binomial distribution with 
"success" parameter eti . 

The method of traffic measurement (on a group) outlined in the pre- 
ceding paragraph has the disadvantage that it collects information very 
slowly. But it is relatively cheap, since all that has to be recorded is 
whether the first trunk is busy at arrival epochs or not, and it has the 
additional advantage that its statistical theory is relatively simple and 
has been well developed in the literature. It must be kept in mind that 
the sampling error estimates we develop are limited to measurements 
made at epochs just preceding arrivals. 

Often the traffic engineer needs to estimate the load offered to a group, 
rather than the load carried by it. The use of S to estimate ai tells him 
what fraction of the time the first trunk is busy. However, there are 
cases in which the knowledge of a x determines the offered load. This 
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Fig. 21 — The covariance value R(l) as a function of traffic a for uniformly 
distributed intcrarrivals. 



occurs, in fact, whenever Oi is a monotone function of the offered load 
a only. For example, when arrivals are Poisson, we have ay = a/{\ + a), 
so it is reasonable to use S/(n — S) as an estimator of the offered load 
a. When arrivals are regular, a,\ = e~ l °, so a reasonable estimate of a 
is l/(log n — log S). 

In the Poisson example, this method of estimating a can be evaluated 
readily if we estimate aT instead by means of (n + l)/(# +1) — 1, 
whose stochastic limit is obviously a~ . 
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Fig. 22 — The covariance value R(2) as a function of traffic a for uniformly 
distributed interarrivals. 
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Fig. 23 — The covariance value i?(3) as a function of traffic a for uniformly 
distributed interarrivals. 
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Fig. 24 — The covariance value R(l) as a function of traffic a for regular 
arrivals. 
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Fig. 25 — The covariance value R(2) as a function of traffic a for regular 
arrivals. 
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Fig. 26 — The covariance value fi(3) as a function of traffic a for regular 
arrivals. 

The generating function of S is 

E{x s \ = [1 + (x - IK]" - Z.r/Pr{£ = j). 



;=(] 



Hence 



E{S + l) -1 } = f E{x s }dx = 
Jo 



E 



S+ 1 



-1 = 



"i 



1 - (1 - Oi) 
(n + l)ai 

[1 - (1 - ai)"]. 



n+l 



There seems to be no simple formula for the second moment of this 
estimator, nor for that of n/S. However, noting that 



[n + 1: 



< 



(n + l) 2 



(£ + l)(5 + 2) - (S + l) 2 ' 
we can verify (by the same method as above) that 

(fi + l) 2 1 n + 1 1 + aid - ai)" + 



/•: 



0S+l)(£+2)j n + 2 



ai z 
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Fig. 27 — R(n) and R(0)ai" as functions of traffic a for Poisson arrivals and 
N = 3 trunks. 



aid so conclude that 



var 



n + 1 
S + 1 



1 ^ 



n + 11 + a,(l - a,)" +1 [1 -(1 -ai)" +1 ] 2 



n + 2 



or 



«r 



This lower bound is likely to be very close to the variance on the left 
for large n, so that, in this region, 

jn + i A ( fll + 2)(i - ai y +1 - (i - ai y-" +2 

var< : — 1 > *~ = . 

l*S + 1 j tt! 2 

It can easily be shown (by the methods of Section IV) that, if N = °o ; 
i.e., if the trunk group is unlimited in size, the covariance function is 
exponential in character: 

R(n) = R(0) ai n . 
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This suggests that, in some cases, R(n) ~ R(0)ai" is a good approxima- 
tion to the covariance for N < °° . This approximation is equivalent to 
ignoring the correction terms Qh in the recurrence relation (12) for the 
covariance. Since the sign of the Qk in (12) is negative, it is clear that 
the approximation is an overestimate. 

The covariance R(n) for n = 0, 1, 2, 3 and the overestimate R(0)a,i" 
for R(n) have been plotted together in Figs. 27, 28 and 29 for 3, 5 and 
8 trunks, respectively, and Poisson arrivals. The curves suggest the fol- 
lowing conclusions: 

i. The approximation R(n) ^ R(0)a" is likely to be good if the load 
per trunk a/N is low. 

ii. If the load per trunk a/N is high, e.g., a/N = 1, the approxima- 
tion R(n) ~ R(0)ai" may give a figure for the covariance (between 



I 1.4 



c_ 1.2 
(0 

o 

en 











R(n) and R(0)aP as 

FUNCTIONS OF TRAFFIC a 
FOR POISSON ARRIVALS 
























r" "*« 














/ 

/ 
/ 
t 

1 ^ 




X 




\R(0) 






1 

/ / 

/ / 


1 / 
/ 

f / 


y" 


^^ 


s 








II 
If 
1 1 

II 


/ 
/ 
1 
1 

1 


s' 


' 


"^N 


•s.n=2 

.J! =3*" 






ll, 


i ^-r 1 

/ / 
/ / 








RU>> 


**>.^ 


i 

/_/ 

7 / 

/ 'I 
1 il 
1 1 1 
I lit 


7 1 

7 /// 


/ 
I 








v?(2) 




w 










s^R(3) 




w 


// 















2 3 4 5 6 

a, OFFERED TRAFFIC IN ERLANGS 



Fig. 28 — R(n) and R{0)a\ n as functions of traffic a for Poisson arrivals and 
N = 5 trunks. 
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Fig. 29 — R(n) and R(0)cii" as functions of traffic a for Poisson arrivals and 
N = 8 trunks. 

separate observations of N{k)) that is several times the actual value. 
This effect seems to increase with the separation, n. 
iii. Variances, such as that of 

S = S N(k), 
i 

computed on the basis of the approximation R(n) ~ R(0)ai" are overesti- 
mates, so that use of this approximation in estimating sampling error is 
conservative. 

iv. The value of a at which R(n) has its (apparently unique) maxi- 
mum seems to be the same for all n, depending only on N, the size of 
the group. 
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